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Epigeal spider responses to fertilization and plant litter: testing biodiversity theory at the ground level 
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Abstract. Recent studies of nutrient additions to terrestrial ecosystems have focused on the “aerial” portion of the food 
web associated with living plants. These studies showed nutrient loading increased arthropod abundance and biomass, but 
decreased diversity. However, none of these studies explicitly examined nutrient loading effects on epigeal arthropods. To 
test nutrient loading effects on epigeal spiders and on individual species within a temperate-latitude grassland community, 
we used pitfall traps to sample spiders for four years within 24 large (314 m 2 ) plots in which wc manipulated nutrients 
(NPK fertilizer) and plant litter (litter removed or left in place). We measured the diversity, abundance, biomass, and 
community structure responses of the spider community, and of wolf spiders (Lycosidae) and linyphiid spiders 
(Linyphiidae), as well as the abundance and biomass responses of the six most common species. We hypothesized increased 
nutrient loading would increase epigeal spider abundance and biomass but decrease diversity. Contrary to predictions, 
spider species richness, diversity, and biomass were not significantly affected by fertilization, while fertilization resulted in 
significantly increased abundance. Also contrary to predictions, plant litter did not affect any of these variables. Linyphiid 
spiders had the strongest responses to fertilization, with significantly increased abundance and biomass, and, contrary to 
predictions, increased species richness in fertilized plots. Wolf spiders responded more closely to predictions. Our results 
indicate that the epigeal spider community does not respond as would be predicted by biodiversity-productivity theory. 

This underscores the need to integrate the largely detritus-based epigeal community into current biodiversity-productivity 
theory. 
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Human activity has resulted in a significant increase in the 
global nitrogen (N) pool through fertilization and increased 
atmospheric N deposition (Vitousek et al. 1997; Fenn et al. 
2003; Galloway et al. 2003). Typical plant community 
responses include decreased plant species richness, increased 
standing crop biomass, and the limitation of community com¬ 
position to a few dominant species (Hector et al. 1999; Tilman 
et ah 2002; Suding et al 2005; Patrick et al. 2008a). This 
research has provided substantive support for biodiversity- 
productivity theory, which predicts declines in local and 
regional richness as one moves from mesotrophic to eutrophic 
systems (Grime 1973; McCann 2000; Worm & Duffy 2003; 
Suding et al 2005; Chalcraft et al. 2008). Similarly with 
terrestrial arthropods, increased nutrient loading has been 
linked with decreased species richness and increased abun¬ 
dance, particularly among those species most closely linked to 
the living-plant portion of the food web (Knops et al. 1999; 
Haddad et al. 2000, 2001; but see Patrick et al. 2008b). This 
“eutrophication effect” (Fenn et a!. 2003) can result in 
significant biodiversity loss and potential decline in important 
ecosystem functions, such as ecosystem stability (McCann 
2000; Larsen et al. 2005). 

Previous studies of nutrient loading have focused on the 
portions of the food web closely tied to living plant material; 
e.g., the “aerial” arthropod community associated with the 
upper portions (e.g., stems) of plants (Knops et al. 1999; 
Haddad et al 2000, 2001). Although there is evidence to 
support the eutrophication effect on aerial arthropod diver¬ 
sity, less is known regarding how nutrient loading affects the 
epigeal (ground-level) arthropod community. A differential 
response by the epigeal arthropod community may result from 
it being more closely tied to the detritus-based portion of the 


food web (Halaj & Wise 2002). Despite the important role it 
may play, the epigeal arthropod community remains an under- 
studied food web component (Wardle 2002; Hattenschwiler 
et al. 2005; Cross et al. 2006). 

Nutrient loading also increases plant litter production 
(Long et al. 2003; Patrick et al. 2008a), which can increase 
the basal food resource for the detrital community and 
increase detritivore and epigeal predator abundances (Halaj 
et al. 2000; Halaj & Wise 2002; Moore et al. 2004). 
Furthermore, plant litter increases habitat complexity, which 
can also increase arthropod abundance and diversity (Lawton 
1983; Strong et al. 1984; Rypstra et al. 1999). Although more 
plant litter production could increase detritivore and epigeal 
predator abundance and biomass (Halaj et al. 2000), a 
reduction in litter diversity could result in decreased diversity 
of detritivores and epigeal predators (Hattenschwiler & Gasser 
2005; Wardle 2006), mirroring the aerial community response 
to nutrient loading. 

Spiders, in particular, are abundant generalist predators 
(Wise 1993) that can significantly impact terrestrial food webs 
(Wise et al 1999), and epigeal spiders (e.g., Lycosidae and 
Linyphiidae) are closely linked to the detritivore community 
(Wise et al. 1999; Chen & Wise 1999; Wise 2006). The 
abundance of epigeal spiders is limited ultimately by the 
abundance of their mainly detritivorous prey via bottom-up 
forces through the detritus-based portion of the food web 
(Chen & Wise 1999; Wise et al. 1999; Wise 2004, 2006). Thus, 
increasing plant detritus can increase spider abundance by 
increasing the quantity of food available to their detritivorous 
prey (Chen & Wise 1999; Wise et al. 1999; Wise 2004). 
Increased detritus also enhances habitat structure for hiding 
and web building (Uetz 1979, 1991; Rypstra et al. 1999), which 
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can also moderately increase the local richness of the spider 
community (Rypstra et al. 1999), also differentially affecting 
individual spider species. 

Even though spider abundance may increase, spider 
diversity may not increase proportionally because the redueed 
diversity of plant detritus can limit predator diversity in the 
detrital food web (Hattensehwiler & Gasser 2005; Wardle 
2006). Thus, it is reasonable to expect that predators 
dependent upon the detritivore food web may have the same 
response to fertilization as predators more closely associated 
with the aerial food web. Even though more plant litter is 
produced, increasing the resource base of the detritivore food 
web, lower litter diversity likely begets lower detritivore and 
detritivore-predator diversity (Hattensehwiler & Gasser 2005; 
Wardle 2006). Interestingly, no epigeal spider studies (focused 
strictly on cursorial spiders; e.g., wolf spiders) have looked at 
the spider diversity response to basal resource manipulation. 
Moreover, no studies have examined responses of the 
predominantly epigeal spider family Linyphiidae (wandering 
sheet/tangle-web builders) that may patrol multiple webs at 
ground level (Uetz et al. 1999). 

Here we report the results of a four-year study that 
investigated the response of the epigeal spider community to 
experimental manipulations of NPK fertilization and plant 
litter availability in a temperate-latitude grassland. We 
measured the diversity, abundance, biomass, and community 
structure responses of the entire epigeal spider community, the 
spider families Lycosidae and Linyphiidae, and the dominant 
individual spider species. Our goal was to integrate the detrital 
food web into biodiversity-productivity theory through insight 
gained from the responses of predators that rely largely on the 
detritivore food web. Based on previous studies that sampled 
the aerial arthropod community responses to nutrient loading 
(e.g., Knops et al. 1999; Haddad et al. 2000, 2001), we tested 
two hypotheses: ( 1 ) fertilization will cause spider biomass and 
abundance to increase and spider species richness to decrease, 
and ( 2 ) the presenee of plant litter will moderately increase 
spider species richness, though this effect will be dampened in 
fertilized plots. 

METHODS 

Study site and experimental design, —The study was done at 
the 163.5 ha Bath Nature Preserve (BNP: 41°10'36.2"N, 
81°38'58.7"W), Bath Township, Summit County, Ohio, USA, 
in a 16 ha section of grassland. Until the early 1980s, the study 
site was a hay meadow, harvested one or many times per year. 
Since then, the area has been mown annually in late August to 
early September, and the mown vegetation has been left on the 
field. The dominant vegetation is an herbaceous, graminoid 
community largely dominated by cool-season C 3 grasses, 
e.g., Bromus inermis Leyss., Lolhim aniuciinaceum (Schreb.) 
Darbysh, Phlewu pratense L., and Anthoxanthum odoratum L. 
The site is moderately productive relative to other grasslands 
within the upper Midwest and across the USA (Patrick et ah 
2008a). The dominant soil type is Ellsworth silt loam (E1B), 
which consists of moderately well drained, moderately deep to 
deep soils formed in silty clay loam or clay loam glacial till of 
the Wisconsin Age (Ritchie & Steiger 1974). 

During August 2001, twenty-four 20-m diameter circular 
plots (314 nr) were established. These experimental plots were 


separated by at least 20 m and were at least 30 m away from 
any other habitat. Treatments were applied in a 2 X 2 factorial 
design of fertilizer (+F = fertilizer added, — F = no fertilizer) 
and plant litter (—L = litter removed, +L = litter left in situ 
after yearly mowing) with control plots characterized as no 
fertilization and plant litter left in situ (+L/-F), resulting in six 
replicates per treatment. Hereafter, all references to “litter” 
refer to the previous year’s mown vegetation and any 
vegetation senesced and found within the sampling quadrat 
after standing crop removal. In April 2002 and continuing 
each April through 2005, Scotts brand Osmocote 8-9 month 
Slow Release Fertilizer 19-6-12 (NPK; Scotts, Marysville, 
Ohio USA) was applied at 20 gNm -2 in fertilized plots, well 
above the Kochy & Wilson (2005) 15 g Nm ~ 2 yr -1 threshold 
necessary to induce a eutrophication effect in grasslands 
and other habitats. We could not exclude ambient wet/dry 
atmospheric N deposition, though deposition rates from 1990 
to 2005 were relatively low at —1.01 g N m “ 2 yr ^ 1 at a nearby 
monitoring site in Lykens (162 km west of our study site), 
Ohio, USA, and —0.93 g N m ~~ 2 yr -1 at another nearby 
monitoring site in Mercer Co. (G.K. Goddard site, 96 km east 
of our study site), Pennsylvania, USA (US EPA 2005). Within 
two days of annual mowing of the whole site by the local 
township with a large tractor and brush hog mower (autumn 
2001-2004), litter was removed from litter-removal treatments 
using a small 23 hp lawn tractor with a pull-behind 8 hp Agri- 
Fab Mow-N-Vac trailer attachment (Agri-Fab, Sullivan, 
Illinois, USA). 

Spider community sampling, —Spiders were collected using 
four pitfall traps in each of the 24 experimental plots (n = 96 
total pitfall traps). Within each plot, a single trap was placed 
5 m from the center of the plot at each of four magnetic 
compass directions (northeast, northwest, southeast, and 
southwest). Each trap consisted of a 10 cm diameter, 18 cm 
tall PVC sleeve into which a 710-mL plastic cup was inserted 
and filled to —4 cm with a 50/50 water/propylene glycol 
mixture. To deter trap raiders (e.g., microtine mammals), 
prevent captured spiders from climbing out of the trap, and 
prevent precipitation from directly Hooding the trap, an 8 -cm 
powder funnel with a base enlarged to — 3 cm was inserted 
and a 15 cm X 15 cm board was placed over each trap, leaving 
— 3 cm clearance. Starting mid to late May (mid-July during 
2004) and continuing through mid to late August, traps were 
alternately left open for two weeks and closed for two weeks. 
This resulted in three sampling periods each year during 2002, 
2003, and 2005. During 2004, only the second and third 
sampling periods were collected. When each two-week 
sampling period ended, the plastic cups were removed, the 
contents collected and preserved in 70% ethyl alcohol, and the 
PVC sleeve was tightly capped. Although pitfall traps do not 
capture all spiders in the community, they are an effective 
sampling technique for determining the relative abundanee 
and species richness of epigeai spiders (Greenslade 1964; 
Phillips & Cobb 2005). Spiders captured in each trap were 
identified to speeies for all mature specimens (taxonomic 
names follow Platniek 2012), and to family for all immature 
specimens, and exact numbers of species/families within in 
each trap were recorded and dried at 70°C for 72 h to 
determine speeies-specific biomasses to the nearest 0.0001 g. 
Lacking sufficient numbers captured within a trap, some 
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extremely small species did not register a biomass, and their 
biomass was recorded as “0.000001 g” to differentiate them 
from true zeroes in analyses. 

Statistical analyses,—We tested the effects of fertilization, 
plant litter, and the interaction of fertilization and litter on the 
abundance, biomass, species richness (SR), and effective 
Shannon’s diversity (e H , where H' is the Shannon diversity 
index) of (1) all mature spiders (Araneae), (2) lycosid and 
linyphiid spiders and (3) abundance and biomass of the six 
most abundant spider species. We used e H to correct for 
differences in species richness that might have resulted from 
differential spider abundances (Ricklefs & Miller 2000; 
Haddad et al. 2000). To calculate the average SR within a 
plot, we summed the total number of spider species caught in 
each trap, then averaged this SR for each of the four traps 
within a plot within a sampling period (including zeroes for 
traps where no spiders were captured), then averaged these 
SRs for each plot across sampling periods in a year, yielding 
n = 24 samples within each year. The same method was used 
to calculate the average abundance, biomass and e H ' within a 
plot within a year, also yielding n = 24 samples within each 
year. Correlations and regressions of these spider responses 
with plant species richness (plant SR) and standing crop 
biomass utilized data from Patrick et al. (2008a). 

To analyze trends per year and per treatment in abundance, 
biomass, SR, and e H ' we used SAS software version 8.01 (SAS 
Institute, Inc. 1999) to calculate a Generalized Linear Mixed- 
effect Model (GLMM) in PROC MIXED with Type III 
effects based upon the covariance structure of compound 
symmetry, and Gaussian distribution of errors. The various 
models used the different response variables (biomass, SR, e H \ 
abundance), and for the predictor variables used fertilized vs. 
unfertilized, litter removed vs. litter left in situ, year, and their 
factorial interactions, with year as the repeated predictor. 
WTien year was detected as a significant effect for a response 
variable, we tested for treatment effects within a year and 
again used SAS to calculate a GLMM in PROC MIXED with 
Type III effects based upon the covariance structure of 
compound symmetry, Gaussian distribution of error, with 
fertilization, litter, and the factorial interaction of fertilization 
and litter as predictor variables. 

To assess treatment effects on aggregate biotic and abiotic 
components in our system, we applied nonmetric multidimen¬ 
sional scaling (NMS: Kruskal 1964) using PC-ORD (McCune 
& Mefford 2006). For 2005, variables used for each of the 24 
plots were average spider species richness per plot and four 
variables used in a previously published analysis (Patrick et al. 
2008a): average plant litter biomass, average PAR per plot, 
average percent soil moisture per plot, and average percent 
soil organic content per plot, resulting in a matrix with five 
columns and 24 rows (plots). The same analysis was run a 
second time with the same variable, except Linyphiidae species 
richness replaced spider species richness. Because (1) NMS is 
scale sensitive, (2) these variables are on radically different 
measurement scales, and (3) variables have an enormous range 
of values between variables, data were transformed to 
proportions relative to the highest value for each variable 
(i.e., each value in a column was divided by the largest value in 
that column, creating a unitless range from 0-1 for each 
column). The NMS analysis was run with Sorensen distance, 


time as the random seed for the starting configuration, 9999 
runs stepping down from 5 to 1 dimensions with the real data, 
999 Monte Carlo runs to assess the probability of a similar 
final stress obtained by chance, and a 0.005 stability criterion. 
Additionally for 2005 and to support NMS analyses with 
stable results, we used PC-ORD (McCune & Mefford 2006) 
to run the multi-response permutation procedure (MRPP: 
Mielke 1984) to test for the hypothesis of no difference among 
treatments. The MRPP used Sorensen distance with the four 
treatments as the a priori groupings, resulting in a matrix with 
five columns (biotic and abiotic variables) and 24 rows (plots) 
and was calculated with all four treatments together, and for 
pairwise comparisons between treatments to test for the 
strength of difference between individual treatments. 

RESULTS 

General trends.—A total of 13,174 spiders from 14 families 
was captured during 14,784 trap nights. Of this total, 2515 
spiders were immature and from 11 families, while the 
remaining 10,659 spiders were mature and from 94 species 
and 12 families (Table 1). Lycosidae was the most commonly 
captured spider family, with 6577 mature specimens (61.7% of 
all mature spiders captured) from 20 species, while Linyphii¬ 
dae was the second most commonly captured spider family 
with 3200 mature specimens (30.0% of all mature spiders 
captured) from 34 species. Together these two families 
represented 9777 (91.7% of all mature spiders captured) 
specimens from 54 species (57.4% of all species captured). 
Spider diversity, corrected for abundance with e H , was 
significantly affected by fertilization and by year, but not by 
litter (Table 2). The factorial interactions between fertilization 
and year, fertilization and litter, litter and year, and the fully 
factorial interaction of fertilization by litter by year were not 
significant (Table 2). Because of the significance of year 
(Fig. 1), we tested for treatment effects on spider diversity 
within each year and by 2005 (Table 3) neither fertilization, 
nor litter nor their interaction was significant. 

Fertilization significantly affected Araneae (all spiders) 
abundance but not Araneae SR or Araneae biomass (Table 2, 
Fig. 2A- C). Moreover, fertilization effects were significant for 
Linyphiidae SR, abundance and biomass (Table 2, Fig. 2D-F), 
as well as for Lycosidae SR and abundance (Table 2, Fig. 2G- 
H), but not for Lycosidae biomass (Table 2, Fig. 21). All response 
variables were significantly affected by year (Table 2), and 
Araneae SR and abundance, and Linyphiidae SR, abundance 
and biomass had significant fertilization and year interactions. 

Neither Araneae SR nor Lycosidae SR were significantly 
correlated with abundance (r = 0.335 and r = -0.190, 
respectively), but Linyphiidae SR was well correlated with 
abundance (r = 0.857; Fig. 3A-C). As with abundance, 
biomass was only correlated in the Linyphiidae SR (r = 
0.629; Fig. 3D-F). Although Araneae SR was not significantly 
correlated with plant SR (r = 0.276; Fig 4A), Linyphiidae SR 
was negatively correlated with plant SR (r = -0.400), and 
Lycosidae SR was positively correlated with plant SR (r = 
0.639; Fig. 4B-C). Araneae SR was also not correlated with 
standing crop biomass (r = 0; Fig. 4D), while Linyphiidae SR 
was positively correlated and Lycosidae SR was negatively 
correlated with standing crop biomass (r - 0.629 and r = 
-0.425, respectively; Fig. 4E-F). 
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Table 1. Total numbers of each family and species of spider captured during the four year manipulative experiment. “+L/—F” represents the 
control treatment of unfertilized plots with litter left in situ, LI— F” represents the unfertilized with litter removed treatment, “+L/+F” 
represents the fertilized with litter left in situ treatment, “ — L/+F” represents the fertilized with litter-removed treatment, and “Total” represents 
the total number caught. Immature spiders were classified as “unidentified,” and family names are in bold. 


Family/Species 

+L/-F 

-L/-F 

+L/+F 

-L/+F 

Total 

Agelenidae 

unidentified 

1 

0 

0 

0 

1 

Araneidae 

unidentified 

1 

2 

0 

0 

3 

Clubionidae 

15 

15 

4 

3 

37 

Chibiomi abboti L. Koch 1866 

0 

0 

1 

0 

1 

Clnbionci hast oil) Gertsch 1941 

6 

9 

1 

1 

17 

unidentified 

9 

6 

2 

2 

19 

Corinnidae 

90 

71 

29 

26 

216 

Castiaiieira gertscbi Kaston 1945 

3 

5 

0 

0 

8 

Castiaueiva longipalpa (Hentz 1847) 

1 

0 

0 

1 

2 

Castiaiieira variata Gertsch 1942 

1 

0 

1 

0 

2 

Mcriola decepta Banks 1895 

10 

13 

3 

6 

32 

Plmirotimpus borealis (Emerton 1911) 

0 

0 

1 

0 

1 

Scotinella britcheri (Petrunkevitch 1910) 

0 

0 

0 

2 

2 

Scotinella fratrella (Gertsch 1935) 

71 

42 

24 

17 

154 

Scotinella wadisonia Levi 1951 

3 

11 

0 

0 

14 

unidentified 

1 

0 

0 

0 

1 

Dictynidae 

Ocarina arcnata Keyserling 1887 

0 

0 

2 

1 

3 

Gnaphosidae 

139 

104 

77 

77 

397 

Drassyllus creolus Chamberlin & Gertsch 1940 

12 

12 

3 

2 

29 

Drassyllus depressus (Emerton 1890) 

63 

57 

20 

30 

170 

Gnapbosa parvula Banks 1896 

43 

19 

44 

37 

143 

Litopyllns temporarily Chamberlin 1922 

0 

0 

2 

1 

3 

unidentified 

21 

16 

8 

7 

52 

Hahniidae 

5 

17 

4 

2 

28 

Neoantistea agilis (Keyserling 1887) 

0 

1 

2 ‘ 

0 

3 

Neoantistea niagna (Keyserling 1887) 

2 

4 

2 

1 

9 

Neoant istea riparia (Keyserling 1887) 

1 

8 

0 

0 

9 

unidentified 

2 

4 

0 

1 

7 

Liny phi idae 

717 

532 

1131 

1058 

3438 

Agyneta sp. 1 

0 

1 

0 

0 

1 

Agyneta sp. 2 

1 

0 

0 

0 

1 

Agyneta sp. 3 

0 

0 

0 

1 

1 

Alloniengea dentisetis (Grube 1861) 

0 

0 

1 

1 

2 

Bathyphantes pallidns (Banks 1892) 

127 

64 

336 

232 

759 

Centromerus cornupalpis (O. P.-Cambridge 1875) 

4 

0 

11 

4 

19 

Ceraticelns similis (Banks 1892) 

0 

1 

1 

0 

2 

Ceratinopsis laticeps Emerton 1882 

0 

0 

0 

1 

1 

Collinsia pluniosa (Emerton 1882) 

31 

20 

150 

192 

393 

Diplostyla concolor (Wider 1834) 

3 

8 

117 

42 

170 

Eridantes erigonoides (Emerton 1882) 

242 

126 

195 

225 

788 

Erigone autnninalis Emerton 1882 

56 

51 

33 

45 

185 

Erigone den tiger a O. P.-Cambridge 1874 

0 

1 

1 

2 

4 

Grammonota gentilis Banks 1898 

0 

0 

0 

1 

1 

Granimonota inornata Emerton 1882 

20 

51 

13 

17 

101 

fslandiana flaveola (Banks 1892) 

11 

12 

8 

7 

38 

Maso sundevalli (Westring 1851) 

1 

0 

0 

0 

1 

Meioneta fabra (Keyserling 1886) 

10 

13 

6 

9 

38 

Meioneta micaria (Emerton 1882) 

5 

4 

2 

0 

11 

Meioneta nniniaculata (Banks 1892) 

85 

52 

61 

72 

270 

Mennessns entoniologicus (Emerton 1911) 

0 

1 

0 

0 

1 

Merniessus jona (Bishop & Crosby 1938) 

9 

9 

8 

2 

28 

Mennessns tridentatus (Emerton 1882) 

3 

2 

1 

3 

9 
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Table 1.—Continued. 


Family/Species 

+L/-F 

-L/-F 

+L/+F 

-L/+F 

Total 

Mennessus trilobatus (Emerton 1882) 

38 

26 

55 

56 

175 

Neriene clatlirata (Sundevall 1830) 

23 

8 

2 

6 

39 

Oedothorax trilobatus (Banks 1896) 

0 

0 

2 

5 

7 

Tennesseellum formica (Emerton 1882) 

0 

1 

0 

0 

1 

Tcnuiplumtes tenuis (Blackwall 1852) 

0 

0 

0 

5 

5 

Walckenaeria directa (O. P.-Cambridge 1874) 

3 

0 

6 

3 

12 

Walckenaeria palustris Millidge 1983 

0 

0 

1 

0 

1 

Walckenaeria sp. 1 

1 

0 

0 

0 

1 

Walckenaeria sp. 2 

0 

0 

1 

0 

1 

Walckenaeria spiralis (Emerton 1882) 

12 

28 

30 

63 

133 

Walckenaeria tibialis (Emerton 1882) 

1 

0 

0 

0 

1 

unidentified 

31 

53 

90 

64 

238 

Liocranidae 

Agvoeca pratensis Emerton 1890 

1 

0 

0 

0 

1 

Lycosidae 

1913 

2021 

2422 

2276 

8632 

Allocosa funerea (Hentz 1844) 

19 

15 

2 

14 

50 

Hogna lielluo (Walckenaer 1837) 

0 

6 

1 

1 

8 

Pardosa milvina (Hentz 1844) 

0 

3 

0 

0 

3 

Pardosa modica (Blackwall 1846) 

2 

1 

6 

3 

12 

Pardosa inoesta Banks 3892 

344 

233 

1444 

1177 

3198 

Pardosa saxatilis (Hentz 1844) 

140 

132 

16 

73 

361 

Pirata aspirans Chamberlin 1904 

0 

1 

0 

0 

1 

Pirata sedentarius Montgomery 1904 

0 

1 

1 

10 

12 

Piratula canadensis (Dondale & Redner 1981) 

0 

1 

0 

0 

1 

Piratula gigantea (Gertsch 1934) 

0 

0 

0 

1 

1 

Piratula insularis (Emerton 1885) 

7 

27 

10 

22 

66 

Piratula minuta (Emerton 1885) 

365 

638 

538 

395 

1936 

Rabidosa punctulata (Hentz 1844) 

2 

3 

0 

0 

5 

Rabidosa rabida (Walckenaer 1837) 

4 

3 

0 

0 

7 

Schizocosa avida (Walckenaer 1837) 

33 

94 

4 

9 

140 

Schizocosa bilineata (Emerton 1885) 

80 

61 

25 

34 

200 

Schizocosa crassipalpata Roewer 1951 

113 

108 

43 

71 

335 

Trochosa ruricola (De Geer 1778) 

24 

16 

36 

37 

113 

Trocliosa terricola Thorell 1856 

18 

29 

42 

36 

125 

Varacosa avara (Keyserling 1877) 

2 

1 

0 

0 

3 

unidentified 

760 

648 

254 

393 

2055 

Philodromidae 

Ebo latithorax Keyserling 1884 

0 

1 

0 

0 

I 

Salticidae 

63 

75 

16 

21 

175 

Ghelna barrowsi (Kaston 1873) 

1 

2 

0 

1 

4 

Ghelna canadensis (Banks 1897) 

4 

3 

0 

2 

9 

Ghelna castanea (Hentz 1846) 

2 

0 

0 

0 

2 

Marpissa lineata (C.L. Koch 1846) 

7 

3 

1 

2 

13 

Mynnaraclme formicaria (De Geer 1778) 

0 

1 

0 

1 

2 

Neon avalonus Gertsch & Ivie 1955 

1 

0 

0 

0 

1 

Neon nelli Peckham & Peckham 1888 

7 

12 

1 

2 

22 

Neon plutonus Gertsch & Ivie 1955 

25 

35 

7 

7 

74 

Sarinda hentzi (Banks 1913) 

0 

1 

0 

0 

1 

Talavera minuta (Banks 1895) 

11 

7 

2 

3 

23 

unidentified 

5 

11 

5 

3 

24 

Tetragnathidae 

51 

46 

16 

22 

135 

Glenognatha foxi (McCook 1894) 

9 

10 

3 

8 

30 

Pachygnatha autunmalis Marx 1884 

17 

14 

9 

6 

46 

Pacliygnatha clerki Sundevall 1823 

0 

0 

0 

1 

1 

Pachygnatha xanthostoma C.L. Koch 1845 

0 

0 

1 

0 

1 

Tetragnatha laboriosa Hentz 1850 

1 

0 

0 

0 

1 

unidentified 

24 

22 

3 

7 

56 
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Table 1.—Continued. 


Family/Species 

+L/-F 

-L/-F 

+L/+F 

-L/+F 

Total 

Thomisidae 

19 

54 

14 

20 

107 

Xysticus bicuspis Keyserling 1887 

1 

0 

0 

0 

1 

Xysticus canadensis Gertsch 1934 

1 

1 

0 

0 

2 

Xysticus fevox (Hentz 1847) 

4 

25 

7 

9 

45 

Xysticus fraternus Banks 1895 

0 

1 

0 

0 

1 

Xysticus luctans (C.L. Koch 1845) 

0 

0 

1 

0 

1 

unidentified 

13 

27 

6 

11 

57 


Spider species-level analyses.—The most abundant spider 
was Pardosa moesta Banks 1892, with 3198 mature specimens 
captured (30.0% of all mature spiders, 48.6% of mature 
lycosids). Pardosa moesta was easily the most captured spider 
in fertilized plots although being virtually absent from 
unfertilized plots (Fig. 5A), and had a strong response to 
fertilization and year, resulting in a significant fertilization by 
year interaction (Table 3). A smaller lycosid, Piratula minuta 
(Emerton 1885), was the second most abundant spider with 
1936 mature specimens captured (18.2% of all mature 
spiders, 29.4% of mature lycosids) but did not seem to 
specifically and consistently respond to a particular treat¬ 
ment (Fig. 5B), though year and the fertilization by litter 
interaction were significant (Table 3). Together, these two 
species accounted for nearly half (48.2%) of all mature 
spiders captured, and over three-quarters (78%) of all mature 
lycosid spiders. 

The third most abundant spider was the linyphiid Eridantes 
erigonoides (Emerton 1882), with 788 mature specimens 
captured (7.4% of all mature spiders, 24.6% of mature 
linyphiid spiders). Similar to Pi. minuta , E. erigonoides did 
not consistently respond to any particular treatment (Fig. 5C), 
but did have a significant response to year in the repeated- 
measures analysis (Table 3). However, the fourth most 


abundant spider, the linyphiid Bathyphantes pallidus (Banks 
1892) with 759 mature specimens captured (7.1% of all mature 
spiders, 23.7% of linyphiid spiders), strongly responded to 
fertilization, year, and the interaction between fertilization and 
year, and responded marginally significantly to litter (Fig. 5D; 
Table 3). Virtually absent from the study site during the first 
two years of the study, B. pallidus became a fairly common 
spider in fertilized plots during the final two years of the study, 
with a weak affiliation to plots where litter was left in situ 
(Fig. 5D). Also virtually absent from the site during the first 
two years of the study, the fifth most abundant spider, the 
linyphiid Collinsia plumosa (Emerton 1882) with 393 speci¬ 
mens captured (3.7% of ail mature spiders, 12.3% of mature 
linyphiid spiders), also strongly responded to fertilization and 
year, which resulted in a significant interaction between 
fertilization and year (Fig. 5E; Table 3). 

One of the most abundant of the larger spiders, the lycosid 
Schizocosa avida (Walckenaer 1837) with 140 mature speci¬ 
mens captured (1.3% of all mature spiders, 2.1% of mature 
lycosid spiders), was virtually absent from fertilized plots 
(Fig. 5F). All variables were significant at the a < 0.05 level 
(Table 3). The species was most commonly captured in litter- 
removed treatments, particularly in unfertilized plots with 
litter removed (—L/—F). 


Table 2.—Results of the GLMM for each response variable. Given are the F values with degrees of freedom and the resulting P values, where 
bolded values indicate significance at ot < 0.0007 (after applying Bonferroni correction). The predictor variables are fertilization (F), litter (L), 
year (Y), and (*) their fully factorial interactions. Species richness is indicated by “SR.” 


Response Variable F L Y F*L F*Y L*Y F*L*Y 
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Figure 1.—Average effective Shannon’s H' (e H ) of spider species 
in each treatment; the letter “a” above a year denotes significance at 
a < 0.05 for fertilization. Open circles (O) and “+L/—F” represent 
the control treatment plots of unfertilized and litter left in situ, open 
triangles (A) and “—L/—F” represent unfertilized and litter removed 
plots, filled circles (#) and “+L/+F” represent fertilized and litter left 
in situ plots, and filled triangles (A) and L/+F” represent fertilized 
and litter removed plots. 

Aggregate ecosystem-level analyses.—The NMS ordination 
with Araneae species richness showed clustering of plots into 
treatments (Fig. 6A), and the ordination axes explained 78.9% 
of the variance, with the first axis explaining 44.8% of the 
variance, and the second axis explaining 34.1% of the 
variance. The final stress = 6.23 with a final instability = 
0.099, and results of the Monte Carlo simulation indicated 
that this stress was less than expected by chance (P = 0.001). 
Following Clarke (1993), a final stress between 5 and 10 was a 
very good ordination and did not present any real risk of 
misinterpretation. The first axis separated fertilized and 
unfertilized plots with high correlations to PAR ( r = 
-0.912) and soil moisture (/• = -0.612) in the direction of 
unfertilized plots, and correlations to species richness (r = 
0.347) and percent soil organic content (/• = 0.446) in the 
direction of fertilized plots, while plant litter biomass was not 
well correlated (r = -0.127). The second axis separated litter 
removed from litter left in situ plots with a strong correlation 
to plant litter biomass (r = 0.916) in the direction of litter left 


in situ plots and PAR (r = -0.529) in the direction of litter 
removed plots, but only weak or no correlations to the other 
three variables: spider species richness r = —0.230, percent soil 
moisture r = —0.144 and percent soil organic content r = 
-0.034. 

Although the NMS ordination with Lycosidae did not 
produce a stable result, the NMS ordination with Linyphiidae 
species richness showed clustering of plots into treatments 
(Fig. 6 R), and the ordination axes explained 92.9% of the 
variance, with the first axis explaining 71.7% of the varianee 
and the second axis explaining 21.2% of the variance. The 
final stress = 8.89 with a final instability = 0.056, and results 
of the Monte Carlo simulation indicated that this stress was 
less than expected by chance (P = 0.002). The first axis 
separated fertilized and unfertilized plots with high correla¬ 
tions to PAR (r = -0.953) and soil moisture (r = -0.491) in 
the direction of unfertilized plots, and correlations to linyphiid 
species richness (r = 0.815) and percent soil organic content 
(r = 0.466) in the direction of fertilized plots, while plant litter 
biomass was not well correlated (r = 0.183). The second axis 
separated litter removed from litter left in situ plots with a 
strong correlation to plant litter biomass (r = 0.942) in the 
direction of litter left in situ plots and weakly correlated to 
linyphiid species richness (/* = —0.359) and percent soil 
organic content (/• = -0.297) in the direction of litter removed 
plots, and no correlations to the other two variables: PAR r = 
-0.039 and percent soil moisture r = —0.034. 

For both NMS ordinations, the separation of plots into 
treatment clusters was supported by MRPP (Table 4). When 
all four treatments were run together, the null hypothesis of no 
difference between treatments was rejected, with high within- 
group agreement and very strong separation between groups. 
Pairwise comparisons of treatments showed that fertilized 
plots, while still significantly distinct, were more similar to 
each other than fertilized treatment plots were to any of the 
unfertilized treatment plots. The same pattern existed for 
unfertilized plots, with strong separation of unfertilized plots, 
yet with lower dissimilarity than when unfertilized plots were 
compared to fertilized plots. As expected, the maximal 
differences occurred when extremes of treatments were paired, 
as in —L/—F vs. +L/+F, and +L/-F vs. -L/+F, indicating that 
“opposite” treatments significantly alter biotie and abiotic 
components of the local habitat. 


Table 3.- Results of the GLMM for six species. Given are the F values with degrees of freedom, and the resulting P values, where bolded 
values indicate significance at a < 0.0012 (after Bonferroni correction). The predictor variables are fertilization (F), litter (L), year (Y), and (*) 
their fully factorial interactions. 


Response Variable F L Y F*L F*Y L*Y F*L*Y 
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G. H. 





Figure 2.—Species richness, abundance, and biomass of all spiders (A-C), linyphiids (D-F), and lycosids (G-I). Definitions of symbols and 
abbreviations for treatments are given in Figure 1, while the letter “a” above a year denotes significance at a < 0.05 for fertilization. 


DISCUSSION 

Spider abundance increased as a result of fertilization, but 
neither spider biomass nor spider species richness was 
significantly affected, and therefore our first hypothesis was 
not supported. This result is contrary to previous studies in 
which arthropod diversity in fertilized plots decreased as 
abundance increased (e.g., Knops et al. 1999; Haddad et al. 
2000, 2001). However, our null result for the overall spider 
community likely resulted from of a canceling effect of the 
responses of the two dominant spider families, the wolf spiders 
(Lycosidae) and the iinyphiid spiders (Linyphiidae). Wolf 
spider abundance was indeed significantly affected by 
fertilization, but wolf spider biomass and species richness 
were not affected. As nearly two-thirds of all spiders captured 
were wolf spiders, the response of this family drove the 
patterns found in the overall spider data, though it should be 


noted that our pitfall trap sampling method may have more 
bias towards wolf spiders due to their cursorial habit. 

On the other hand, fertilization increased the abundance, 
biomass, and species richness of Iinyphiid spiders during the 
final two years of the study. The species richness response was 
completely opposite of the predictions of biodiversity-produc¬ 
tivity theory (Haddad et al. 2000, 2001; Suding et al. 2005). 
This was likely a result of a bottom-up food web response to 
fertilization by the main Iinyphiid spider food source, 
collembolans (Harwood et al. 2001; Romero & Harwood 
2010). Collembolans did not respond to the treatments despite 
the increased abundance of plant litter (L. B. Patrick, 
unpublished data), with the basal resource for the collembolan 
prey being the bacteria that aid in the breakdown of plant 
litter. Thus, while the basal resource likely increased, the 
primary consumer of that resource, collembolans, did not, but 
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Figure 3.—Regressions of all spiders, linyphiids, and lycosids (left to right) against abundance (A C) and biomass (D-F). Symbols are defined 
in Figure 1. and data presented are for 2005. 
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Figure 4. Regressions of all spiders, linyphiids, and lycosids (left to right) against plant species richness (A-C) and standing crop biomass 
(D F). Symbols are defined in Figure 1, and data presented are for 2005. 


the primary collembolan predator did increase in abundance and 
diversity. It is therefore feasible to propose that a top-down food 
web effect by linyphiid spiders limited collembolan abundance, 
ultimately enhancing their own abundance and diversity. 

The differences in responses between wolf and linyphiid 
spiders are likely the results of different foraging behaviors. 


Spider guilds are based primarily on foraging behavior, 
habitat preferences and web type (Uetz et al. 1999). Based 
upon this classification system, wolf spiders are considered 
ground-running spiders, and linyphiid spiders are character¬ 
ized as wandering sheet/tangle web weavers. Wolf spiders are 
classic epigeal wandering spiders that actively hunt for prey 
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A. Pardosa moesta B. Piratula minuta 




C. Eridantes erigonoides D. Bathyphantes pallidus 




E. Collinsia plumosa F. Schizocosa avida 




Figure 5.—Average abundance of selected species by year. Definitions of symbols and abbreviations for treatments are given in Figure 1, 
while the letters above each year denote significance at a < 0.05 for “a” = fertilization, tk b” = litter, and “c” = the interaction of fertilization 
and litter. 
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Figure 6.—Two-dimensional ordination of ecosystem-level properties from 2005 from nonmetric multidimensional scaling (NMS) using plant 
litter biomass, PAR, percent soil moisture, percent soil organic content, and (A) Araneae species richness or (B) Linyphiidae species richness 
(Lycosidae did not produce a stable result). Vectors indicate the direction and strength of correlations between axis scores and emergent 
properties (/?“ cutoff for joint biplot = 0.100), and ordinations are rotated to the dominant axis of fertilization. The percent of variance explained 
by each axis is noted next to the axis title. See Figure 1 for key to treatment symbols. 
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Table 4.—Results of Multi-response Permutation Procedure 
(MRPP) on emergent properties for 2005 to support NMS analyses 
(Fig. 5). T describes the separation between groups (dissimilarity), 
and A is the chance-corrected within-group agreement. “All” 
indicates all four treatments included in the MRPP, and the 
remainders are MRPP pairwise comparisons of treatments to assess 
dissimilarity (lower T and higher A). +L indicates litter left in situ , — L 
indicates litter removed, +F indicates fertilization, — F indicates 
no fertilization. 


Groups 

T 

A 

P 

All Spiders 

All 

-9.892 

0.421 

< 0.0001 

+L/-F vs. -L/-F 

-3.352 

0.164 

0.0097 

+L/-F vs. +L/+F 

-5.838 

0.332 

0.0006 

+L/-F vs. -L/+F 

-6.014 

0.414 

0.0010 

-L/-F vs. +L/+F 

-6.470 

0.424 

0.0006 

-L/-F vs. -L/+F 

-5.643 

0.313 

0.0009 

+L/+F vs. -L/+F 

-4.999 

0.313 

0.0023 

Linyphiidae 

All 

-10.055 
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< 0.0001 

+L/-F vs. -L/-F 

-3.236 

0.159 

0.0102 

+L/-F vs. +L/+F 

-6.395 

0.373 

0.0005 

+17-F vs. -L/+F 

-6.445 

0.442 

0.0008 

-L/-F vs. +L/+F 

-6.501 

0.435 

0.0006 

-L/-F vs. -L/+F 

-5.842 

0.338 

0.0010 

+L/+F vs. -L/+F 

-5.188 
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and are largely restricted to hunting in two-dimensional space. 
Thus, habitat structure (i.e., the physical structure of the 
surrounding environment, including plant litter and living 
plant material) may enhance wolf spider hunting success by 
providing additional hiding places for ambush hunting and for 
lairs (Rypstra et al. 3999; Halaj et al. 2000). However, too 
much habitat structure can also increase predation and 
intraguild predation risk through increased density responses 
to habitat structure while inhibiting movement (Wise 2006; 
Rypstra et al. 2007). These factors reduce both abundance and 
hunting success, thereby reducing the numbers of wolf spiders 
in an area. In our study, the increased habitat structure that 
resulted from fertilization seemed to moderately (but not 
significantly) reduce wolf spider species richness while 
increasing wolf spider abundance. However, this increased 
abundance was likely due to the population explosion in 
fertilized plots of the medium-sized wolf spider, P. moesta (see 
below). 

Linyphiid spiders rely more upon webs for prey capture, 
sometimes maintaining and patrolling multiple webs (Uetz 
et al. 1999). Although these webs are generally constructed at 
or close to the ground level, the webs can enhance prey capture 
space to include a portion of a third dimension. Moreover, 
increased habitat structure can provide additional structure 
for web building (Rypstra et al. 1999). Thus, while fertilized 
plots significantly reduced wolf spider species richness, 
probably due to the enhanced habitat structure that impeded 
foraging ability, these plots may have provided the tiny web 
building linyphiid spiders the habitat structure to flourish 
because of the increased structure for web building and, thus, 
increased prey capture rates. 

A medium-sized wolf spider, P. moesta thrived in fertilized 
plots probably due to decreased intraguild predation by larger 


wolf spiders (e.g., S. avida) that became less abundant in 
fertilized plots. Moreover, there was increased abundance of 
potential prey in fertilized plots (Patrick et al. 2008b). These 
two factors together likely released P. moesta from competi¬ 
tion and predation, resulting in an increased abundance in 
fertilized plots. However, these microhabitat changes likely 
caused the decreased abundance of S. avida in fertilized plots, 
as the increased habitat structure likely impeded this species’ 
foraging abilities. Piratula minuta was one of the smallest wolf 
spiders captured at our site, and could have benefited from the 
increased habitat structure in a similar way to P. moesta. 
Although it was the second most abundant spider captured 
during the course of our study, it did not significantly respond 
to a single habitat, except for unfertilized plots with litter 
removed during 2005. Interestingly, Pi. minuta was observed 
on several occasions in the jaws of P. moesta and S. avida , 
making Pi. minuta a victim of intraguild predation. 

Spider species richness was not significantly correlated with 
either plant species riehness or standing crop biomass. 
However, both dominant spider families responded to 
fertilization (by the fourth year of the study) in distinctly 
different ways. Wolf spiders followed predictions of current 
biodiversity-productivity theory, with decreased species rich¬ 
ness associated with decreased plant species richness and 
increased standing crop biomass. Although fertilization 
increased wolf spider abundance, wolf spider species richness 
was correlated with plant species richness and therefore 
decreased as nutrient loading into the system increased. 

Finally, despite documented effects of increased habitat 
structure on arthropod abundances and diversity (e.g., 
Lawton 1983; Halaj et al. 2000), particularly for spiders 
(e.g., Uetz 1991; Rypstra et al. 1999; Halaj et al. 2000), our 
results do not support our second hypothesis. Plant litter had 
no significant effect on spider species richness. Most studies of 
spider responses to plant litter have been conducted in plant 
monocultures in agroecosystems (e.g., Rypstra et al. 1999). 
These managed ecosystems tend to have much higher 
disturbance and more bare ground than would be expected 
from a grassland. Thus, increased refugia via plant litter 
additions to these agroecosystems would certainly provide 
more habitat than the existing bare ground, so it is perhaps 
not surprising that there have been stronger responses to plant 
litter in agroecosystems. 

Analysis of the spider community and associated abiotic 
variables demonstrated strong treatment effects. These highly 
differentiated treatments are likely to have a strong effect on 
ecosystem properties (e.g., nutrient cycling, carbon sequester¬ 
ing), an effect likely to increase through time as the treatment 
plots further mature. Spiders have been shown to affect 
detritivore abundance (Wise et al. 1999), thereby indirectly 
altering nutrient cycling within the system (Chen & Wise 
1999). The results of our ordinations clearly showed that our 
plots responded to our treatments and that the spider 
community affected ecosystem-level processes. The long-term 
implications are unknown, but it is clear that the trajectories 
of each treatment are significantly different and may impact 
ecosystem function and services. To our knowledge, this is the 
first time that these biotic and abiotic factors have been 
coupled in a multivariate ordination to explicitly determine 
whether they can define discrete and distinct predator 
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communities and their associated abiotic properties in the 
context of the biodiversity-productivity theory. Most previous 
work (e.g., Haddad et al. 2000, 2001) did not attempt to 
associate the invertebrate community with abiotic changes 
resulting from fertilization, and we know of no other studies 
that coupled fertilization and plant litter effects to test 
predictions of biodiversity-productivity theory. 

The diversity and community structure of spiders and other 
arthropods are sensitive to plot size (Martinko et al. 2006). 
The large size of our experimental plots integrated important 
determinants of the within-plot plant communities, including 
spatial heterogeneity (De Boeck et al. 2006), leaching of 
nutrients from litter (Berendse 1998), local nutrient cycling 
(Hooper & Vitousek 1998) and the translocation of nutrients 
within clumping and clonal plants (Hutchings & Bradbury 
1986), which are the primary growth forms of our dominant 
graminoids (Patrick et al. 2008a). These spatial factors are also 
important to epigeal spiders because of their vagility and their 
need to find suitable food; the larger plot sizes more 
realistically emulate natural habitat patches of varying quality 
and can support higher insect diversity (Martinko et al. 2006). 
Other studies that examined the effects of nutrient loading on 
arthropod communities had plot sizes ranging from 9 m 2 - 
169 nr (e.g., Knops et al. 1999; Haddad et al. 2001), making 
our experimental plots (314 m 2 ) nearly twice as large—an 
important factor when considering the vagility of some spider 
species. 

However, we realize that our study has some distinct 
differences when compared to previous work. Our use of an 
NPK fertilizer, as opposed to N-only fertilizer, is likely to have 
induced a stronger response to fertilization due to the added P 
and K. Nevertheless, our plant results (see Patrick et al. 2008a) 
were generally consistent with other plant studies that used 
NPK fertilizers (e.g., Carson & Barrett 1988; Turkington et al. 
2002) and N-only fertilizers (e.g., Haddad et al. 2000; Tilman 
et al. 2002), which allowed us to formulate our epigeal spider 
hypotheses on the same bases as previous studies that 
investigated the responses of arthropods to nutrient loading. 
Further, our running definition of litter (see Methods) includes 
the vegetation mown in the previous year and not removed 
from litter left in situ treatment plots, potentially altering the 
nutritional quality of the litter relative to naturally senesced 
vegetation, and the physical structure of the litter as it lay 
after mowing (e.g., Semmartin et al. 2004). Because the timing 
of mowing was determined by the local township, litter 
from the annual mowing accumulated earlier than might 
normally be expected for this region of the USA. However, 
were the mowing to stop, the site would very quickly yield 
to encroaching woody vegetation typical of early secondary 
succession. 

Our study underscores the disjunct between conventional, 
plant-based biodiversity-productivity theory and the animal 
component of the food web, particularly epigeal predators. 
This portion of the food web is more closely associated with 
the quality of its basal resource (plant litter) than with the 
diversity of that resource (Cross et al. 2006; Seeber et al. 2008). 
This starkly contrasts with the more aerial portion of the food 
web that is more dependent on living plants, where specialist 
herbivores can be affected by plant diversity more than by 
plant quality. Ultimately, the loss of plant species with 


increased nutrient loading may result in the loss of arthropod 
herbivores and their specialist predators and parasites. 
However, the increases in diversity may be balanced by the 
epigeal community and its different resource base. 
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